Experimental Tomographic State Reconstruction of Itinerant Microwave Photons 
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A wide range of experiments studying microwave photons localized in superconducting cavities 
have made important contributions to our understanding of the quantum properties of radiation. 
Propagating microwave photons, however, have so far been studied much less intensely Here we 
present measurements in which we reconstruct the Wigner function of itinerant single photon Fock 
states and their superposition with the vacuum using linear amplifiers and quadrature amplitude 
detectors. We have developed efficient methods to separate the detected single photon signal from 
the noise added by the amplifier by analyzing the moments of the measured amplitude distribution 
up to 4th order. This work is expected to enable studies of propagating microwaves in the context 
of linear quantum optics. 



The quantum properties of microwave frequency pho- 
tons localized in space are typically investigated in the 
context of cavity quantum electrodynamics (QED) exper- 
iments with Rydberg atoms [TJ [2] or superconducting cir- 
cuits [3H5]. Propagating (itinerant) individual microwave 
photons have been generated [6] but their quantum prop- 
erties have not been studied with the same intensity yet. 
This is partly due to the fact that efficient single photon 
counters in this frequency domain are still under devel- 
opment [7]. However, recently it was shown that char- 
acteristic quantum properties of propagating microwave 
photons, such as antibunching, can be observed in corre- 
lation measurements using linear amplifiers and quadra- 
ture amplitude detectors [EJ [9]. In the context of cir- 
cuit QED [3j [4] , propagating microwaves are also used to 
control [10] and read out the quantum state of artificial 
atoms [IT] and to observe phenomena such as resonance 
fluorescence [12]. 

In general, the quantum state of any field mode a is 
characterized by its density matrix p or an equivalent 
quasi-probability distribution such as the Wigner, the 
Husimi Q or the Glauber-Sudarshan P function [13] [14]. 
Less widely appreciated, the mode a is also equivalently 
specified by the infinite set of its moments ((a^) n a m ) [T5] . 
In this work we use measurements of such field moments 
up to 4th order to characterize single photon states. 

In the optical domain many experimental techniques 
exist to reconstruct the quantum state of light using 
square law detectors [16]. In many instances, homodyne 
detection schemes are used to record the statistical prop- 
erties of single quadrature components for different local 
oscillator phases which allows to reconstruct the Wigner 
function by an inverse Radon transformation [17] . Alter- 
natively, in heterodyne detection schemes the joint statis- 
tics of two conjugate quadrature components described 
by the Husimi Q distribution are measured [18]. Both 
techniques thus allow for the full reconstruction of the 
quantum state of a single field mode. 

In the microwave domain linear amplifiers are used to 
measure the amplitude of the signal instead of its in- 
tensity. In this case, one can make use of phase sensitive 



amplifiers in homodyne detection, ideally amplifying only 
one quadrature of the signal noiselessly, or of phase in- 
sensitive amplifiers in heterodyne detection, amplifying 
both quadratures of the signal equally while adding at 
least the vacuum noise to the signal [19] [20] . 

Here we demonstrate photon state tomography using 
a phase insensitive amplifier in combination with hetero- 
dyne detection of both field quadratures. We first discuss 
the relation of the single propagating field mode to the 
resonator mode acting as the photon source (see Fig. [T]) . 
Then we describe a method to systematically separate 
the relevant signal from the noise using a single linear 
amplifier chain. Similar methods have recently been dis- 
cussed for a setup with two linear amplifier channels [21] . 
Finally, we use this method to reconstruct the Wigner 




FIG. 1: Simplified schematic of the experimental setup. An 
optical analog of the photon source realized in our experi- 
ments is shown as a single sided cavity with one highly reflec- 
tive and one partially transmitting mirror. A single photon 
is generated in the cavity by preparing the qubit in its ex- 
cited state when detuned from the cavity and then tuning it 
into resonance with the cavity for half a vacuum Rabi period 
27r/g. The photon is then emitted at a rate n into the out- 
put mode awt resulting in a exponentially decaying envelope 
of the single photon pulse while the mode a- in remains in the 
vacuum state. The signal is then amplified with effective gain 
G and noise in mode h is added. The amplified signal is down 
converted in a microwave quadrature mixer using a local os- 
cillator (LO). The two quadrature amplitudes X and P are 
recorded using a analog-to-digital converter and stored in real 
time in a two dimensional histogram using field programmable 
gate array (FPGA) electronics. 
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function of a single photon Fock state and its superposi- 
tion with the vacuum state. 

We realize a single photon source in a circuit QED 
setup employing a transmission line resonator of fre- 
quency v r ~ 6.77 GHz coupled to a single transmon qubit 
with vacuum Rabi rate 2g/2n « 146 MHz [9]. By tuning 
the qubit prepared in a superposition state a\g) + f3\e) 
into the resonator for exactly half a vacuum Rabi pe- 
riod we generate single photon states a\0) [3 [22]. 
The state preparation time is short compared to the cav- 
ity decay time r = 1/n « 40 ns. We repeat the photon 
generation every 800 ns which allows us to prepare ap- 
proximately 4 x 10 9 single photon states per hour. 

The field generated in the resonator then decays into 
the output mode a out related to the resonator mode 
A by the input-output boundary condition a out (t) = 
y/RA(t) —a,' m (t) [23] where ai n is in the vacuum state, see 
Fig. [T] We integrate the output signal over a weighted 
time window f(t) to define a single time independent 
mode a = J dtf(t)a out (t). Considering the resonator dy- 
namics A(t) = e~^A(0) + y^e - ^ dre^ 1 a- m (r) [24] 
the choice f(t) = v^e - ^0(t), where Q(t) is the Heavi- 
side step function, leads to the identity a = A(0). There- 
fore a is in the same quantum state as the resonator at 
the preparation time t = 0. 

To characterize the state of mode a we first pass the 
signal through a phase insensitive amplifier chain with ef- 
fective gain G which introduces an additional noise mode 
h. The amplified signal is then split equally into two 
parts and mixed with an in-phase and out-of-phase local 
oscillator, respectively, to simultaneously detect the con- 
jugate quadrature components X and P. Consequently, 
X and P are related to a and h by 

VG(X + iP) = VG(a + tf) = S (1) 

where we define the complex amplitude operator S [8]. 

If the amplifier is quantum limited, i.e. the noise mode 
h is in the vacuum state, one can show that the quantum 
state at the amplifier output can be expressed in terms of 
the Husimi Q function Q out (y/Ga) = ^Qm(d) [25]. This 
implies that the Q function remains invariant under the 
amplification process up to a scaling factor. Therefore 
the measurement results of the operator S are distributed 

as Qout{S). 

The best commercially available amplifiers for frequen- 
cies below 10 GHz have a noise temperature of about 
^noise = 2 K and therefore are far from being quantum 
limited. In this case, the noise mode h is in good ap- 
proximation in a thermal state which is represented by a 
Gaussian phase space distribution. As shown in Ref. [26] , 
the measured distribution 

dM {s) = hJ A2 p p a(P)Qh{s*/^G - n (2) 

of S at the amplifier output can then be interpreted as 



the convolution of the P function of mode a and the Q 
function of the noise mode h. 

We store the results of repeated measurements of S 
in a two-dimensional histogram with 1024 x 1024 entries 
which corresponds to a discretized version of the proba- 
bility distribution D^ P \S). To extract the properties of 
mode a alone we perform both a measurement in which 
mode a is left in the vacuum serving as a reference mea- 
surement for the noise where P a (0) = S^(p) is a two- 
dimensional Dirac J-function resulting in the distribution 

M 

£>[I°X°I](S) = ±Q h (S*/VG), (3) 

and a measurement in which the state of interest |^), 
such as a Fock state |1), is prepared. The measured his- 
tograms for both the vacuum Z)[I ^ I] and for a Fock state 
jjil^iM] are dominated by the noise added by the ampli- 
fier, see Fig.[2|a) and (b). However, when calculating the 
numerical difference of both histograms, see Fig.|2|c), we 
already clearly observe the circular symmetric character 
of the single photon phase space distribution. The small 
deviation from an ideal circular symmetry is explained 
by a slight coherent admixture of the vacuum |0) to the 
single photon Fock state |1) due to small errors in the 
state preparation. 

To further analyze the data we calculate the moments 

^) n S m ) p = J d 2 S (S*) n S m DW(S) (4) 

of the two histograms up to a given order that is chosen 
to be n + m = 4 in our experiments. This specific choice 
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FIG. 2: (a) Measured quadrature histogram D mm (S) for a 
in the vacuum where S — \[G (X + iP) . The inset shows a 
horizontal cut through the histogram (red). The distribution 
is well described by a normal distribution (blue) with width 
a = 5.7 (indicated by black arrows) corresponding to a system 
noise temperature of T no ise ~ 21K. (b) Quadrature histogram 
D^^\S) for preparation of single photon Fock states, (c) 
Difference of the two histograms L> [|1)(1|] (5) and D m{ ° l] (S). 
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will be justified later in the text. When noise and signal 
are uncorrected the calculated moments correspond to 
the operator averages 

{(&) n S m ) p = 

G ^ S (^) (7)<(at)V><^(fct r -i>, ( 5 ) 

reducing to ((&) n S m ) m{ol = G^ {h n (h^) m ) when a is 
in the vacuum state. Equation ([5| can then be inverted to 
calculate the moments ((a^) n a m ) from ((S^) n S m ) p and 
((S'"'") n 6' m )|o)(o| up to the desired order, 4 in this case, 
as shown in Fig. |3ja). We note that the quadrature 
histograms are normalized such that the zeroth order 
moments are always unity for all prepared states. The 
off diagonal elements in the moment matrix express co- 
herences between different photon number states. They 
vanish for states with circular symmetric phase space dis- 
tributions such as pure Fock states or thermal states. For 
the Fock state |1) (Fig. J3^a)), we observe that all off di- 
agonal moments are close to zero. In addition, we note 
that the fourth order moment ((a^) 2 a 2 ) is also close to 
indicating antibunching of the prepared single photon 
states [9 j. In contrast, a thermal state with the same 
mean photon number would also have vanishing off di- 
agonal moments but finite diagonal 4th order moments. 
Experimentally, for the single photon Fock state, the 
aforementioned residual coherent admixture of the vac- 
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FIG. 3: Absolute value of the normally ordered moments 
\{(a^) n a m ) \ up to 4th order for (a) a single photon Fock state, 
(b) a superposition state (|0) — |l))/>/2, and two coherent 
states with amplitude (c) a = 1, and (d) a = 0.5. 



uum state leads to a non- vanishing small mean amplitude 
| (a) | = 0.044 and a slightly reduced mean photon num- 
ber (a^a) = 0.91. For an integration time of 12 hours for 
each state, we find errors of the 4th order moments to be 
approximately ±0.1 where the statistical error in the mo- 
ments is known to increase exponentially with increasing 
order [8]. In comparison, the estimated statistical er- 
rors for the first, second and third order moments are 
approximately 1.5 x 10 -3 , 4.5 x 10 -3 , and 1.5 x 10 -2 , 
respectively. 

We have also prepared and analyzed superposition 
states of the type (|0) + e^|l))/>/2, see Fig. ^b). The 
relative phase <p is controlled by the phase of the corre- 
sponding qubit state that is mapped into the resonator. 
For this class of states, the mean amplitude ideally equals 
the mean photon number \(a}\ = (a^a) = 0.5. The first 
equality remains valid even if the state is slightly mixed 
with the vacuum. We have been able to use this property 
to determine the effective gain G of our amplifier chain 
because first and second order moments have a different 
characteristic scaling with G. This allowed us to scale 
X and P axes of the histograms (Fig. [2| such that they 
correspond exactly to real and imaginary part of a-\-hf. 
From our measurement data, we extract \(a)\ = 0.466 
which is close to the expected value. 




FIG. 4: Wigner function W(a = X + iP) for (a) a single 
photon Fock state and (b) a superposition state both recon- 
structed from the measured moments shown in Fig. [3] 
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To further confirm the validity of our scheme, we have 
generated coherent states \a) with amplitude a = 1 and 
a = 0.5 by applying 10 ns square coherent pulses with 
controlled amplitude to the weakly coupled input port 
of the resonator. The moments of coherent states are 
given by ((a^) n a m ) = (a*) n a m . For a = 1 all moments 
are observed to be close to 1 (Fig. |3jc)), as expected. 
This also demonstrates that systematic errors in the de- 
tection chain, such as small nonlinear it ies, are negligible 
as all moments take their expected values. For a = 0.5 
(Fig. |3jd)), the measured moments decay exponentially 
with ((a+) n a m ) = 0.5 n+m , as expected. The fourth order 
moments appear larger than the third order ones, due to 
their larger statistical error. 

From the measured moments we have reconstructed 
the Wigner function W(a) for a single photon Fock state 
and its superposition with the vacuum (Fig. [4|. It is 
sufficient to evaluate 1 



W(a) =J2f 



d 2 A ((gl) n a m )(\*) n \ m e _ 
7T 2 n\m\ 



i|A| 2 +aA*-a*A 



(6) 



up to order n + m = 3 because ((a^) 2 a 2 ) ~ 0. In gen- 
eral, all higher order moments with n + m > 2N — 1 
have to be zero if one diagonal moment ((a^) N a N ) van- 
ishes, which follows from the fact that diagonal moments 
(k\(aJ) n a n \k) with n > k are zero for Fock states \k). 

The Wigner function of the single photon Fock state 
Fig. |4^a) shows clear negative values which indicate the 
quantum character of the observed state. The slight shift 
of | (a) | = 0.044 from the phase space origin that we al- 
ready observed in the raw measurement data (Fig. [2^) of 
the |1) state is also apparent in the reconstructed Wigner 
function. The superposition state ( 1 0) — 1 1 ) ) / y/2 displayed 
in FigQb) has a finite mean amplitude which leads to 
the finite center of mass of the distribution. Still, neg- 
ative values in the distribution persist, illustrating the 
quantum coherence between the |0) and |1) state. We 
have also varied the relative phase <\> of the superposition 
states and have observed the according rotation of the 
Wigner function (data not shown). 

In summary, we have reconstructed the Wigner func- 
tion of itinerant single microwave photons and small co- 
herent fields using linear amplification, quadrature am- 
plitude detection and efficient data analysis even in the 
presence of noise added by the amplifier. We have imple- 
mented a method to separate the quantum signal from 
the amplifier noise in a measurement setup with only 
one detection channel. We believe that propagating mi- 
crowaves will be investigated more intensely in the con- 



text of future quantum optics and also quantum informa- 
tion processing experiments [27] where low noise para- 
metric amplifiers [28H30] have the potential to signifi- 
cantly improve the detection efficiency. 
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